The Josephson current in Luttinger liquid-superconductor junctions 
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We study the Josephson current through a Luttinger liquid in contact with two superconductors. 
We show that it can be deduced from the Casimir energy in a two-boundary version of the sine- 
Gordon model. We develop a new thermodynamic Bethe Ansatz, which, combined with a subtle 
analytic continuation procedure, allows us to calculate this energy in closed form, and obtain the 
complete current-crossover function from the case of complete normal to complete Andreev reflection. 



Low-dimensional condensed matter systems exhibit an 
exciting set of unusual properties, which have been ac- 
tively studied theoretically for some time, and have re- 
cently become a topic of explosive experimental interest. 
The physics of these systems is crucially related to in- 
teractions, which in low dimensions behave nonpertur- 
batively, giving rise to such nonintuitive phenomena as 
spin and charge separation, or charge fractionalization. 

Dealing with these interactions theoretically is a chal- 
lenge, especially as far as transport properties are con- 
cerned. These are of course the most relevant from an 
experimental point of view. New methods have had to 
be developed over the years to tackle these issues, like 
bosonization, conformal field theory and integrability. 
The area where most progress has been made is probably 
that of single impurity problems, including the Kondo ef- 
fect edge states tunneling |3j, and quantum dots. 

One of the paradigms of low-dimensional electronic 
systems is the Luttinger liquid, and it can only be ex- 
pected that fascinating properties should appear when 
combined with another great paradigm of solid state 
physics: superconductivity. Indeed, it was recently 
shown that, as a result of the low-dimensional interac- 
tions, junctions between superconductors and Luttinger 
liquids behaved very differently [|]-[7| from the usual ones 
in particular, there is now an RG flow as the temper- 
ature or the length of the junction is changed, leading, 
for repulsive (resp. attractive) bulk interactions to per- 
fect normal (resp. Andreev) reflection at low energies 0. 
In the non-interacting case there is no flow, and the rel- 
ative amounts of normal and Andreev reflections are set 
by parameters like the superconducting gaps only. 

The formalism developed so far allows one to predict 
global features of the RG flow, but not to exactly com- 
pute physical quantities like the current. This is due to 
the presence of Luttinger interactions, which require a 
non perturbative approach. In this letter, we present a 
new formalism to deal with this difficulty, and show how 
it gives rise to closed form expressions for the current, all 
the way from the UV to the IR fixed points. 

We use here the formulation developed in Q. Fig- 
ure 1 illustrates our considerations. A quantum wire is 
suspended between two bulk superconductors Si r . We 



describe Si and S r with BCS theory, using gaps A\ 



(bulk) 



and A r bulk) e l 



with A 



(bulk) 
Lr 



real numbers. The phase 



difference between the gaps, Xi is crucial in driving the 
Josephson current. Since the electronic states in Si and 
S r are suppressed by the bulk BCS gaps, an effective 
theory at low energies (compared to the gaps) can be 
obtained by integrating out all states in Si and S r . 
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FIG. 1. The effective theory for a Josephson junction with 
a quantum wire between two superconductors Si and S2 , with 
bulk BCS gaps A^J, : a bounded Luttinger liquid, with 
BCS-like couplings A;. r living on the boundaries only. 



The only contributions left from the original supercon- 
ductors are two boundary BCS terms modifying the Lut- 
tinger liquid Hamiltonian: in terms of fermionic left- and 
right-movers, we end up with the boundary contributions 

H B oc A,*UoW (0) + Are*** (R) - 



: HTV-"-y*ij.v-"-; -T- — (1) 

in which ... means h.c. + L <-> R. The values of the 
boundary parameters A/. r are related to the original bulk 

gaps a[ 6 " ZA ^ , and the physical processes at the bound- 
aries depend on them. Namely, two things can happen 
when an electron hits one of the boundaries: it can be 
normal-reflected as an electron, or Andrccv-rcflectcd as 
the time-conjugate excitation (a hole), thereby creating 
an additional Cooper pair in the condensate. The rel- 
ative amounts of normal versus Andreev reflections at 
each boundary is set by A: for A = 0, only normal reflec- 
tion occurs, while for A — > 00, only Andreev reflection 
occurs. The Josephson current can thus be physically 
understood as excitations being Andreev-reflected back 
and forth between the two boundaries, thereby transfer- 
ring charge from one end to the other. 
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Upon bosonization, the theory becomes H = Hq + Hb 
where Hq is the charge part of the Luttinger liquid Hamil- 
tonian, and the boundary Hamiltonian involves the dual 
charge boson (since we don't consider any nontrivial 
spin-scattering processes at the boundaries, the spin and 
charge sectors decouple. We thus consider the charge sec- 
tor only, and drop the spin boson in all further formulas) : 

H B oc Aj cos |£.(0) + A r cos (~4> C (R) - x). (2) 

In the above notation, the interaction parameter is re- 
lated to standard bosonization conventions according to 
f3 2 = (in) 2 R 2 = A/g p = 2/K p . The current through 
the system is then easy to calculate, being given by the 
difference between right- and left-mover densities, which 
under bosonization becomes I(x) — — e J2a(^Ra^ — 
^ La ^Lcr) —ed x (f) c . Reabsorbing this linear derivative 
in the quadratic form of the Hamiltonian then shows that 
the zero-temperature current is given by the simple rela- 
tion I(x) — 2eJf-E (Ai, A r , %), where E is the ground- 
state energy of the system. 

We show in this letter how the Josephson current can 
be obtained exactly using the effective two-boundary 
model described above. The first step, generalizing the 
proposal of |?J , is to think of the bulk theory as the mass- 
less limit of a sine-Gordon theory for the dual charge bo- 
son <j) c , and to use the fact that the boundary sine-Gordon 
terms in (||) do not destroy the integrability of the sine- 
Gordon model ||. Thus, we can derive our results from 
a study of the massless limit of the double-boundary sine- 
Gordon model 



S 



,11 i dx\ [(d^) 2 - (d x 



m cos / 



dt 



At cos^(0)+A r cosf|^)-x 



(3) 



As proposed in || , the key role of integrability is to pro- 
vide one with a basis of the free boson Hilbert space 
made of excitations that scatter in a simple, factorized 
way, with no particle production, both in the bulk and at 
the boundaries. In contrast with ||] and most other work 
on boundary integrable field theories, we have however to 
deal with a problem involving two boundaries. This gives 
rise to unexpected subtleties. 

In imaginary time, two quantization schemes are pos- 
sible 0. In the first one, called the L— channel, we 
consider an inverse temperature L — > oo and write the 
partition function as 



Tr e 



-LH, 



(4) 



where the trace is over eigenstates of the Hamiltonian 
with appropriate boundary conditions at the left and 
right ends. Alternately, we can permute the definitions 
of space and imaginary time and move to the R— channel 



picture, where time flows from one "boundary" state to 
the other. In this case, the partition function is given by 

,LH \Br) (5) 



Z=(B,\e 



where the Hamiltonian is quantized along length L with 
simple periodic boundary conditions. Thus, the states 
\Bi tT ) become initial and final states, linked through time 
evolution along a slice of length R. The advantage of this 
approach is that the interpretation of the boundary states 
is physically clear: as no momentum can flow through 
the boundaries, only states of zero total momentum can 
be emitted or absorbed. As the boundary sine-Gordon 
is an integrable theory |J, factorizability of the scatter- 
ing imposes very strong constraints on the shape of the 
boundary states, constraints which are in fact sufficient 
to determine them completely. 

We specialize in this paper to the free fermion point 
1 = 47r, which corresponds to R 2 C — 1/47T, in the bulk 
attractive regime of the Luttinger liquid. All difficulties 
due to the presence of two boundaries appear in this case, 
while the bulk scattering is trivial and does not obscure 
the issues. There are two fundamental excitations in the 
bulk: the soliton and the antisoliton, carrying rapidities 
9 (their creation operators being denoted A\ 2 (9)), with 
bulk scattering matrix S = — 1. The boundary scattering 
however is nontrivial, and involves nondiagonal reflection 
processes. This affects the boundary states, whose gen- 
eral form is that of a coherent state of excitations with 
opposite rapidities, i.e. || 



\Bi ir >cx exp 



d0K$(O)Al(-O)Al(O) |0 



> 



(6) 



with amplitudes Kf h r (9) explicitly given in JT^]. 

Given the two boundary states, the partition function 
can be calculated in the following way. For a complete 
set of states {a}, we have 



(a\a) 



(7) 



Now it is clear that the only states having nonzero inter- 
nal product with the boundary states are those involving 
pairs of excitations: if rapidity 9 is occupied, so must 
rapidity —9. As we are going to do thermodynamics us- 
ing these states as a basis, it is important to classify 
them properly. Using only positive rapidities for the la- 
beling, we can construct four different pairs, which we 
label 11, 22, 12, 21. A pair 12 has for example the mean- 
ing of a soliton at rapidity —9 and an antisoliton at ra- 
pidity 9. A difficulty then occurs when one implements 
fermionic statistics. Whereas fermions of the same type 
cannot occupy the same rapidity state, we can in fact 
superimpose some pairs, e.g. pairs 11 and 22 but not 
12 and 22. To do the thermodynamics, it is useful to 
consider that in fact none of these pairs can be superim- 
posed, and then account for the possible superpositions 
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by introducing one additional species, a quartet, obtained 
by superposition of 11 and 22, or, equivalently, 12 and 
21. This appears in the boundary state with amplitude 
KftKf* - K} 2 K^ r = det(K hr ). The trace over states a 
now becomes a trace over configurations of these five ex- 
citations on the lattice of allowed rapidities. The fugacity 
of a as determined by the boundary states reads 

(Bi\a)(a\B r ) 



(a\a) 



N 



il[det[Ki(0])K r (9])} 

i=i 



(8) 



where in this example we take a to have N pairs of types 
eij, bi at rapidities &i, i = 1, N, and N q quartets at ra- 
pidities 6f,i = l,...,N g . The TBA can then be done Q 
according to standard procedures by looking for a saddle 
point in the sum (^) , yielding an expression for In Z which 
we can reinterpret as the zero-temperature ground-state 
energy when L — > oo. The end result is 



= lim 

m^O 27T 



d6hx 



l + tr[KiK r ]e 



-2mR cosh 



+ 



det[KiK r ]e 



-4mR cosh 9 



icoshfl. (9) 



Taking the limit, and substituting for the K matrices 
from |12j] . this formula becomes 



in fact in the variables \hi\, \h r \, and careful manipula- 
tions allow one to reproduce ( |l2| ) in the case hih r > 0. 
Expression ( pi] ) cannot be right however since it does 
not distinguish (even after the non-analyticity in hf r has 
been taken into account) the cases hih r > and hih r < 0. 
These cases should be different, and we think ( pi] ) applies 
to the former only. A quick assessment of the situation 
can be made in the limit of very large boundary fields, 
where h\h r > corresponds to H — h fixed boundary con- 
ditions, where based on conformal field theory we expect 
and do indeed get from ( |TT| ) = -=^, while the other 
case describes H — fixed boundary conditions, for which 
we expect Eq = + ^ , the difference being the bound- 
ary dimension of the spin operator, A = 1/2. 

The correct expression for hih r < is obtained us- 
ing the following physical considerations. The point is 
that, for a given boundary field h, the boundary state 
is fully defined only if one specifies the boundary condi- 
tions asymptotically far away from it. The results in || 
assume implicitly that the spin at infinity is fixed at + 
while the boundary field is positive. The case of a spin 
fixed — at infinity should be obtained by considering an 
'excited' boundary state. Following this procedure here 
leads to the result that, for hih r > 0, 



E%(h,h')-EZ(h,- 



-hf) = — k+ 
' 2 + 



(13) 



dn In 



AAfAj cos2x 



-kR 



2A?)(«- 



2A? 



(K-2Af)( K -2A2) _ 2Kl 
(k + 2A2)( k + 2A2) 



(10) 



In fact, expression ( pj| ) hides unexpected subtleties. To 
understand the latter, and how to deal with them, it is 
useful to discuss for a while a simpler and closely related 
model, the critical Ising model with boundary magnetic 
fields hi,h r . As discussed in ||, this theory is integrable, 
and the same kind of thermodynamical approach gives 
rise to the ground-state energy in a finite size R pj: 



x In 



1 



E I (h,h') = - 
-8wh 2 R e 



1 



%Tl{h') 2 R 



de x 



e + 8nh 2 R e + 8i:(h') 2 R 



(11) 



This is in apparent contradiction with physical expecta- 
tions. Indeed, it involves only the squares of the magnetic 
fields, while basic symmetry considerations lead one to 
expect an energy that depends on the fields hi , h r them- 
selves; lowest order perturbation theory in fact leads to 

El{h, h ') = ^- f( h2R ) - f( h ' 2R ) - 2M ' + - (I 2 ) 

with f{x) = 2x 2 (cst. + hix)/R. The contradiction is 
partly resolved by realizing that expression (pi]) expands 



where k + is a particular solution of the quantization con- 



dition 



a iRk 



S-Kh^+ik 
8irh' 2 -ik 



%ir(h'Y+ik 
S-Kih'Y^-ik 



= — 1. Specifically, all so- 



lutions are analytic functions of hf and h 2 , except for 
one pair going as k± ~ ±87r|/iz/i r | for small fields. It can 
be proven that equation ( Jl3|) does indeed provide the an- 
alytical continuation of (O) into the region hih r < 0. 
For instance, the lowest term in the perturbative expan- 
sion reads now for hih r < : — 2n\hih r \ + iir\hih r \ = 



2tt \hih r 



-2irhih r , as desired. We also observe that at 



large magnetic fields, k + — > ^, producing the expected 
gap in the conformal limit. Note that a similar analytical 
continuation procedure was proposed in Jl3[ 

Expression (^) is extraordinarily simple; while each 
of the terms in the left-hand side requires calculation 
of a TBA integral, their difference is simply obtained 
by the solution of an elementary quantization condi- 
tion. We have checked its validity by computing numeri- 
cally the ground state energies of the critical Ising model 
with boundary magnetic fields, and studying their scal- 
ing limit. Results for hi — ±h r will be presented in figure 
^, and are in very good agreement with our formulas. We 
note that the lattice Ising model with boundary magnetic 
fields has never been solved to the best of our knowledge 
(the continuum version with hi — h r has been treated in 
fuj ) , except in the case hi — |H| , where one can check 
that the results agree with those of our TBA approach. 
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FIG. 2. Ground-state energy difference for the 2-boundary 
Ising model with hi — ±h r , for system sizes N = 5 to N — 10, 
and compared to the theoretical prediction (solid line). 

In the sine-Gordon case, we have a similar situation. 
Since Eq in ( JlCj ) is an even function of x, we can con- 
centrate on the region \ > 0. The crucial observation 
is that ( |i~0| ) is valid only in the domain \ € [0, 7r/2[ 
since it exhibits a singularity at x = tt/2, where the 
argument of the log has a zero that hits the real axis 
at k = 0. Better intuition is gained by studying the 
limits A^ r — > oo, where the integral reduces to i?o = 
— irS Jo°° ^ m [l + 2 cos 2xe _:E + e~ 2a: ] . This integral is 

2 

tabulated, and we get Eq = — ^fjj + for x £ [0; 7r /2], 



and E = 
J cx 



MR t 2^ for X S [7r/2,7r]. The current 
^ thus experiences a discontinuity at x — 7r /2, 
which contradicts perturbation theory: it should be odd 
in x, and a smooth function over the interval [0, 7r]. As 
in the Ising case, our procedure to repair this starts 
by identifying the root ko of the quantization condi- 



tion 2 



[(ikf+gfgl cos 2 X ] 



-2ikR 



[ik+gf][ik+gl] ~ C + [*+5f][ife+ff; 

which requires analytic continuation beyond x = tt/2. 
The continuation of the 2-boundary sine-Gordon ground- 
state energy to all values of x then reads 



[ifc— gi][ifc— g;] „-4ifc.B 



1 



E^(A l ,A r ,x)=E , 
E SG (A l ,A r , x )=E -k , 



\X\ < tt/2, 
tt/2 < Ixl < 7T. (14) 



The limiting behaviours of the boundary contributions to 
the ground-state energy are SEq g — > for A/, r large, 
and SEq G oc — A;A r cosx for small A/ jr . 

With this formula, we can write down the Joscphson 
current of our device at the free fermion point 1 = 47r, 
in terms of the dimensionless parameters A\ R = RAf r : 



T , . 8eA? A„ . „ 
Hx) = ^sin2 X 

JO 



dn x 



2A|] 



2A|] 



^+4A|A|cos2 X (k-2A|)(k-2A|) 
1 ^ Z (k+2A^)(k+2A|) L ^ ( K +2AJ)( K +2A^) e 



(15) 



with the understanding that this should be analytically 
continued beyond x = 7r /2 using the above procedure. 
This formula for the current provides complete interpo- 
lation between the limits of the sawtooth function for 
perfect Andreev boundary scattering (A/. r — > 00), and 
the sinx behaviour for small Andreev scattering. These 
and further results are discussed in more detail in 
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FIG. 3. The Josephson current (in units of R 1 ) as a func- 
tion of x, for RA = 0.5, 1, 2, 4, 8 and 16 (in ascending order). 
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